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Abstract 


An  extremal  principle  for  accounting  balance  of  a resource  value  - 
transfer  economy  of  W.  P.  Drews  is  developed.  Through  *t . ex'Stence 
and  uniqueness  are  completely  characterized  in  terms  of  an  associated 
linear  programming  problem  of  pure  distribution  (transportation)  type. 

An  existence  theorem  of  Dantzig's  (established  by  the  Brouwer  fixed 
point  theorem  under  economically  ouestionable  conditions)  is  an  immediate 
special  corollary.  The  principle  is  equivalent  to  unconstrained  minimi* 
/ation  of  a simple  strictly  convex  function  and  can  be  computed  thereby. 
Alternately,  it  is  an  extended  geometric  programming  problem  whose 
dual  is  minimization  of  a strictly  convex  function  subject  to  the  same 


distribution  constraints  mentioned  above. 


1.  Introduction 


Events  of  recent  years  in  re  exhaustible  resources  led  W.  P.  Drews 

to  formulate  a model  fwhich  we  have  not  seen)  in  which  an  economy  is  driven 

by  a transfer  of  agreed  monetary  values  of  r resources,  e.  g.  oil,  food, 

labor,  etc.,  into  agreed  monetary  receipts  of  s re  source -owner  ( - good  s - 

consumer)  groups  such  that  (1)  total  resource  valuation  equals  total  monies 

received.  (2)  prices  of  resources  and  sizes  of  owner  groups  adjust  to  drive 

consumer  goods  prices  and  industrial  activity  levels  ' employment  levels) 

so  that  the  value  of  each  resource  equals  the  sum  the  economy  spends  for 

it  and  the  receipts  of  each  owner-group  equals  the  sum  of  its  expenditures. 

The  latter  detail  (2)  we  call  "accounting  balance.  " 

A variant  of  Drews'  model  was  presented  and  studied  by  Dantzig  in 

[1|.  He  shows  only,  by  making  a strong  (and  economically  questionable) 

assumption  ona  matrix  relating  resource  use  to  owner  group  iize 

and  employing  the  Brouwer  fixed  point  theorem,  that  relative  prices  and 

relative  group  sizes  exist  which  satisfy  the  accounting  balance.  No  possi- 

* 

bilities  of  computation  are  suggested,  the  question  of  uniqueness  is 
untouched,  and  nothing  is  mentioned  concerning  the  effect  of  the 
structure  of  the  matrix  or  the  agreed  monetary  divisions  on  the  balance. 

In  the  following  we  completely  characterize  such  an  economy  by  a 
non-linear  extremal  principle  which  may  be  interpreted  as  the  minimization 
of  an  economic  potential  function.  Accounting  balance  occurs  if  and  only 
if  there  exists  a point  of  economic  equilibrium  (minimum  potential).  Despite 


.‘J . 


the  non-linearity  of  the  functional  minimized,  the  effects  of  the  K*y  (rx  s) 

matrix  \1,  the  resource  valuations  vj vr  and  the  owner  receipts  6j, . . , 6 

onth*  nature  of  the  minimum  (non-existence,  existence,  uniqueness)  are  com - 
pletelv  characterized  in  terms  of  an  associated  linear  programming  prohl*  m of 

distribution  (transportation)  type  with  supplies  Vj vr  receipts  6S 

and  with  the  non -zero  entries  of  M designating  the  possible  "transfer''  routes. 
The  extended  dual  theorem,  p.  182, |2|,  at  1 dual  considerations  are  the  key  here. 

The  extremal  problem  is  re. lured  to  the  unconstrained  minimization 
of  a strictlv  convex  function  involving  onlv  exponential  and  linear  terms. 
Computation  is  easily  possible  with  SHMT,  Piecewise  linear1  zat'on,  or  anv 
other  standard  method.  Alternately  't  can  be  viewed  as  an  extended  geometric 
programming  problem,  whose  dual  is  a convex  programming  problem  with 
pure  distribution  constraints.  The  latter  formulation  seems  most  advantageous 
to  compute  large  systems. 

2.  Notation  and  Formulation 

l.et  ,A,  Aj  be  respectively  the  ith  row,  jth  column  of  the  matrix  A. 

T 

!.et  p , q be  the  row,  column  vector  of  prices  per  unit  resource,  sizes  of 
owner  (consumer)  groups;  y.x  the  vectors  of  goods  prices,  industry  activity 
(employment)  levels.  R,  the  (rxs)  matrix  of  resources  per  unit  industrial 
activity.  C the  (mxs)  matrix  of  goods  per  unit  owner  (consumer)  group 
size,  N the  (nx  m)  matrix  of  industrial  activity  per  unit  good  (in  ( 1]  the 
inverse  of  a Leontief  type  matrix);  \1  * If  \C.  We  assume  Vj,  flj  > 0 lor 
all  i,  j.  Note  that  It,  N,  C are  non-negative  matrices. 

(2.  1)  y r*  pTRN  . x*  NCq 

(2.2)  . v = ' 6 i.e.  the  total  value  is  transferred 
' J J 

The  accounting  balance  conditions  are, 

(2.3)  Vj1  Pj(jMq),  6j-(p*Mj)qj  , all’i.j 


j 
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l.cmma  1 If  v^,  fij  > 0.  for  all  i.jifidCL3)  holds,  then  p.  seui,  qj*evJ  > 0,  for 
all  i.  i,  and  M has  a positive  entr>  in  each  row  and  each  column. 
Lemma  2~  At  accounting  balance, 


v - p ^\1q  = 6 

1 1 J J 


Then  (3)  may  be  rewritten 


(2.4)  ''  M eU«  ♦ VJ  = v 

J ,J 


T.  Mj.eU*  * VJ  - 6, 
i J J 


. all  i.  j 


Let 


(2.5) 


(u,  v)  * T,  M 


«.  J 


'J 


u.  ♦ V 
e i vj 


vu 

II 


V* 

J 


Then,  by  the  Kuhn-Tucker  Theorem,  since  £(u,  v)  is  convex  and  analytic 
in  (u.  v). 


Theorem  1 ^(u.  v)  has  a minimum  at  (u.v) 

if  and  only  if 

the  accounting  balance  (2.4)  holds 


The  problem  of  balance  may  therefore  be  rendered  by  the 
Kxtremal  Principle 

(2.5)  minimize  p1  Mq  - T.  Yjlnpj  - T.  6^  Inq^ 


T 

with  p , q > 0tor  equivalently  * 

u,  ♦ v. 


(2.6)  minimize  5(u,  v)  * £ M e ' VJ  - T V:U:  -T.  6 v 
x i.j  'J  i 1 1 j J J 

without  constraints  . 

Note:  (2.5)  may  be  recognized  as  an  extended  geometric  programming  problem 


3.  Existence 


The  necessary  condition  on  M in  Leirma  1 is  not  sufficient  even 
that  C (u.  v)  be  bounded  below.  To  characterize  this  situation,  consider  the 
following  "transfer"  system  (of  "pure  distribution"  type  in  linear  programming 
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jarpon) 

x v,  i -I.  . . . , r 
je.Hi)  'J  ' 

(3.  1)  v . . 

.r!(j)X'J  “j*  J* 8 

. i 0, 

'J 

where  Mi)  * {j  M >0},  Mj)  £ {»:  M-  > 0} 
rheorem  2 £ (u,  v)  is  bounded  below 

if  and  only  if 

the  transfer  system  (3.  1)  has  a solution 


Proof  If  (3.  1)  has  a solution  x-^,  we  can  rewrite 

(3.  2)  E (u,  v)  - v M , 'VI  xnlu'v) 
r «.J  ‘-1  i,  j M ' J 

and  Ret 

(3.3)  C(u.v>*  S x, , f 1 - ln(x, , / M 1 1 

( i,  j )«P  ,J  ,J  ,J 

where  P * f <i . j ) * x > 0} 

•J 

If  (3.  1)  has  no  solution,  we  minimize  £ 0*  x;,  subject  to  the  transfer 

•.J  'J 

conditions.  The  dual  theorem  { 2]  implies  that  a sequence  (un.  vn)  exists 
for  which  ^ (un.  vn)  -•  . 

% 

Q.  K.l». 


I.et  I)  he  the  coefficient  matrix  of  the  left  side  of  the  equations  in 
(3.  1).  Its  maximum  possible  rank  is  r*s*l.  By  the  regularization"  pro- 
cedure of  Charnes  and  Cooper  |2|  all  cases  can  be  reduced  (also  computational! v ' ' 
to  this  case.  We  discuss  this  elsewhere. 

Here  we  note  that  the  existence  of  a solution  to  (3.  1),  while  sufficient 
for  an  infimum,  is  not  sufficient  for  existence  of  a minimum  of  C(u,  v). 


* 


-5- 


Theorem  3 £<u,  v)  has  a minimum  (-accounting  balance) 

if  and  only  if 

the  transfer  system  (3. 1)  has  a solution  x^  whose  non-zero  values 
designate  a sub-matrix  of  I)  of  equal  rank  to  IJ 


Proof  Suppose  (3.  1)  has  a solution  whose  non-zero  components  designate 

a submatrix  of  rank  less  than  1).  Ix*t  w.  =u.  *v  for  x . > 0.  Let 

•J  ' J 'J 

Zkl*  Uk * V1  ^kl  ® and  *'Mn  tf,e  are  independent  variables 

i.e.  the  uk  * Vj  cannot  be  expressed  (linearly)  in  terms  of  the  w(j. 

We  can  now  write 

H.4)  (?(u,v)  1 C(w.z)  * T.  (M.  c '■*  - x..w  ) ♦ T.  M.,e  ** 

r I„  >0  'J  'J  'I  lk.ll/  K kl 


where  K * {(k.1)  Mk|  >0,  xk|"0}. 

Clearly,  ^(w.  7)  < ^ (w.  z)  whenever  z * z and  some  IM<*W  Th“* 

g(u.v)  has  an  infimum  and  no  minimum  whenever  K i $.  If  K-#,  then 

(3.5)  £(u.v>*  J (MiieW|  1-x.jW.  ) k r •(Mllewij-x  w.  ) « £(w) 

Xjj  > o J ,J  ,J  *.J)€Q  J 'J  'J 

where  x.  ^ satisfies  (3.  1 ),  Q * {fi.  )>•  x.  >0]c{(i,j):x>0}  and 

*•  * 'J  ^ 

A 

Q designates  linearly  independent  columns  of  D,  hence  independent 
variables  w^. 

Rut  £<w>  as  |w|  Since  if  and  only  if  the  rank  of  the  submatrix 

of  1)  designated  by  x^  > 0 equals  that  of  I),  the  theorem  follows. 

Q.  K.  I). 


Dantzig's  result  is  the  corollary. 
Corollary  1:  If  M > 0,  ^(u,  v)  has  a minimum*. 
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C oroll.try  2 If  I’  * ( <t . j>  If  > 0)  • Q • {(i,  j):  M(J  >0}  and  the  columns 
of  D designated  by  P are  linearly  independent,  then  the  minimum 
of  f»(u,  v)  takes  place  at  w^  = Uj+Vj-  Infxjj/Mjj),  (i.j)  e P. 


•4.  Uniqueness 

T 

Assuming  a minimum  exists,  if  p ,q  is  a (price,  size)  solution  so 
■r  T _ - 1 

is  p r ip  , q - (a  >q  for  any  n > 0.  I'niqueness  is  possible  at  most  to 


rotative  prices  and  sizes.  However,  extending  the  argument  in  Theorem 


3.  we  can  show  there  exists  a "basic"  solution  Xjj  of  (3.  1)  so  that 


C.i.  6) 


xn  'J 


x .>  0 
•J 


where  f^iw)  s linear  in  the  independent  variables  w(^.  Evidently 
the  right  hand  is  a strictly  convex  function  of  the  w,j.  (Without  assuming 
a minimum.  £(u,  v)  may  be  written  in  form  (3.6)  plus  possible  addtional 
exponential  terms  in  i.e.  as  a strictiv  convex  function  of  w^  and  zkj. ) 
Theorem  4:  The  minimum  of  £(u,v)  is  unique  in  variables  Wjj  of  (3.6). 
Theorem  5:  If  w(J  minimizes  (3.6),  then  the  set  of  minimizing  (u,  v)  is 
the  set  of  all  solutions  to 


(3.  7) 


u.  ♦ 


V 


•J 


> 0 


u ♦ v s f (w  ) , x = 0 
‘ J »J  »J  »J 


where  the  submatrix  of  D associated  with  tlie  $,s  is  a basis  for  (3.  1). 

Questions  of  uniqueness  and  stability  of  the  system  to  arbitrary  but 
small  perturbations  of  the  v.,  6^  are  closely  related.  Stability  of  solution 
to  such  perturbations  implies  that  the  rank  of  D must  be  r+s-1,  e.  g. 
Theorem  6 If  no  subsums  of  the  and  6j  are  equal,  then  for  a minimum  to 
exist  the  rank  of  I)  must  he  r*s-l.  The  minimum  point  >s  unique  up 
to  relative  prices  and  sizes. 


J 
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Sometimes  M may  be  decomposable  i.e.  block  diagonal  on  suitable 
interchanges  of  rows  and  of  columns,  ff  there  are  B blocks,  the  problem 
splits  into  B independent  problems  of  the  same  form.  When  soluble,  each 
independent  problem  minimum  can  be  at  most  unique  up  to  a scale  factor 
for  prices  (and  its  reciprocal  for  sizes).  It  must  be  so  unique  if  the  block 
problem  is  stable  under  local  perturbations  of  its  own  v|(  6^. 

5.  Computation 

The  accounting  balance  prices  ami  sizes  and  corresponding  consumer 
goods  prices  and  employment  levels  may  he  computed  from  (2.  6)  and  its 
variants  by  any  of  the  standard  convex  programming  methods,  e.  g.  SUMT 
or  piecewise  linearization.  Alternately,  Newton  - Ha  phson  schemes  can  be 
applied  to  solve  the  balance  equations  directly,  now  that  one  knows  the 
theory  of  their  ambient  system.  A third  dire.  non.  which  may  be  most 
convenient  for  large  systems  is  through  solving  the  geometric  programming 
dual  to  (2. 5). 

From  (3 1,  page  210  et  jeq.  the  extended  geometric  program 

9 % 

T *V  1 J/ - ei  t •&«  T 

<5.  1)  minimize  p Mq  + ln(pj  . . . Prrq,  ‘-  .qs  ) with  p . q > 0. 

has  a dual 

(5.2)  minimize  -v(6)  * T 6 fln(6j./e)  - InM  | 

(|.J)CQ,J  ,J  'J 

subject  to  (3.1)  on  the  6tJ. 

If  6^  is  an  optimal  solution,  then 

(5.3)  6^  « M(jp*qj  . (i.j)  e Q. 

Note  how  (5.3)  recovers  the  result  in  (3.71. 
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Since  linear  programming  problems  of  distribution  type  may  be 

computed  at  least  100  times  faster  than  general  linear  programs,  e.  g. 

2 

in  a few  seconds  on  existing  codes  for  r.  h - 10  (and  larger),  and  since 
existing  codes  for  distribution  problems  handle  by  far  the  largest  linear 
programs,  this  mode  of  solution  seems  most  attractive  for  large  systems. 

6.  Extensions 

The  model  needs  extension  in  the  directions  of  more  explicit 
connections  between  employment  levels,  consumer ‘group  sizes,  etc. 

Multi -p«Tiod  extensions  would  also  be  desirable,  as  would  introduction  of 
stochastic  (risk  or  uncertainty)  elements.  The  authors  are  currently 
making  such  extensions. 


*' 
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